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Abstract. We study a family of differential operators L a in two 
variables, depending on the coupling parameter a > that appears 
only in the boundary conditions. Our main concern is the spectral 
properties of L a , which turn out to be quite different for a < 1 and 
for a > 1. In particular, L Q has a unique self-adjoint realization for 
a < 1 and many such realizations for a > 1. In the more difficult 
case a > 1 an analysis of non-elliptic pscudodiffcrential operators 
in dimension one is involved. 



1. Introduction 

In the paper [8] Smilansky suggested a mathematical model which 
he called "The irreversible quantum graph". In this model a one- 
dimensional quantum graph interacts with a finite system of harmonic 
oscillators attached at different points of the graph. Regardless of the 
physical meaning of this model, it is quite interesting from the math- 
ematical point of view, since, being a singular perturbation problem, 
it exhibits many unusual effects. These effects appear already in the 
one-oscillator case. They were discussed in the survey paper [9], see 
also references therein. 

In the simplest case (the graph is a real line, with only one oscillator 
attached) the problem consists in the study of a family of differential 
operators A a on M 2 , depending on the coupling parameter a > 0. The 
differential expression which defines the action of A a does not involve 
a, this parameter appears only in the transmission condition across 
a straight line in the plane. The operator A admits an exhaustive 
description via the separation of variables, and the passage to A a with 
a^O can be expressed, at least formally, in the terms of perturbations 
of quadratic forms. The main peculiarity of the problem stems from the 
fact that the perturbation is too strong: it is only relatively bounded 
but not relatively compact with respect to the operator Ao (in the sense 
of quadratic forms). For this reason, the standard machinery of the 
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perturbation theory does not work. Still, it turned out to be possible 
to give a detailed description of the spectrum <t(A q ) for all a > 0. A 
borderline value a* of the parameter a exists, such that the properties 
of o~(A a ) are quite different for a < a* and for a > a*. For a < a* 
the absolutely continuous (a.c.) spectrum of A a is the same as for 
A , including the multiplicity. Eigenvalues appear below the bottom 
of er(Ao), their number grows indefinitely as a S a* and satisfies 
an asymptotic relation of a non-standard type. For a = a* these 
eigenvalues disappear and a new branch of the a.c. spectrum appears 
instead, filling [0, oo). For a above the threshold a*, the operator A a 
is not semi-bounded any more and its a.c. spectrum fills the whole 
real line. Thus, the system exhibits a sort of phase transition as the 
parameter a crosses the threshold a*. 

The mathematical mechanism behind such a behaviour of the spec- 
trum lies in a very special form of the transmission condition for the 
operator A a . This condition generates in a natural way an infinite 
Jacobi matrix which depends on the parameter a and whose spectral 
properties for a < a* and for a > a* are quite different. 

The papers [4] and [5] are devoted to the case of two oscillators, 
but actually their results show what happens in the general case of an 
arbitrary number of oscillators. It was an initiative of Des Evans, to 
start the work on these papers, and we take pleasure in emphasizing 
his role in the study of this class of problems. 

In the present paper we investigate another family of differential op- 
erators, say L a , of a similar nature. It was also proposed by Smilansky 
(private communication). Again, all operators in the family are de- 
termined by a differential expression not depending on the parameter, 
and they differ by the transmission condition. Like in the case of the 
family A a , a certain family of Jacobi matrices is closely related to the 
operator. However, the properties of the two families are rather dif- 
ferent and another type of phase transition occurs. Namely, for large 
values of a the operator L a has many self-adjoint realizations, and the 
negative spectrum of each realization is discrete and unbounded from 
below. The mechanism of this transition lies in an unusual breaking of 
the Shapiro - Lopatinsky ellipticity condition in several points on the 
interface line, and the analysis of this situation involves a study of a 
priori estimates for some non-elliptic pseudodifferential operators. 

In the last section of the paper we briefly consider yet another family 
M a of differential operators. It looks rather similar to the family L Q , 
but some important details in the behaviour of the spectrum are quite 
different. 
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Taken together, the families A a , L a , and M a show that presence 
of the coupling parameter in the boundary condition may cause quite 
different types of the phase transition. It is tempting to develop a 
general scheme which would include all these examples as special cases. 

2. Stating the problem. Preliminaries 

We study a family L a of differential operators on the cylinder Q = 
Kx §' identified with the strip R x (0, 2n) with periodic boundary 
conditions for all functions involved. Further on, x stands for the co- 
ordinate on R and y for the co-ordinate on S 1 . The operator L a is 
generated by the Laplacian —AU = —U" 2 — U' y 2 and two conditions at 
x — 0. The first condition is the continuity 

(2.1) U(0+,y) = U(0-y) (=U(0,y)) 

and the second one is a 'transmission condition' at x = 0: 

(2.2) U' x (0+, y) - U' x (0- y) = ia {U' y (0, y) cosy + (U(0, y) cosy);) . 

In (2.2) a is a real parameter. The passage a i— > —a corresponds to 
the change of variables y i— > y + tt, which does not affect the spectrum. 
For this reason it is enough to consider a > 0. 
By using the Fourier expansion 

(2.3) U=(2n)-^J2 u ^ x y ny 

(in short, U ~ {u n }), we reduce the problem formally to an infinite 
system of ordinary differential operators on the real axis, 

(2.4) -AU ~ {- u » + n 2 u n }, x^O, n e Z, 
coupled by the conditions 

(2.5) u n (0+) = u n (0-) (= « n (0)) , 

(2.6) <(0+) - <(0-) = -a((n + l/2)« n+1 (0) + (n - l/2)u n ^(0)). 

The operator L is just the standard Laplacian on the cylinder Q, 
with the domain H 2 (Q). Thus, for a = the above formal reduction 
of the partial differential operator is legal, the system decouples, and 
we get 

Here — J^- stands for the self-adjoint operator in L 2 (R) with the domain 
H 2 (R) and the symbol J^ ffi denotes the orthogonal sum of operators. 
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The expansion (2.7) leads to the complete description of the spectrum 
<t(L ): it is absolutely continuous, fills the half-line [0, oo), and its 
multiplicity function is given by 

(2.8) m a . c .(A;L ) = 2 + 4[A], VA > 0, 

where, as usual, [A] denotes the integer part of a real number A. 

For a ^ we must first specify in what sense the conditions (2.1) 
and (2.2) are understood. Suppose that U G L 2 (Q) is a weak solution 
of the equation —AU = F G L 2 in each semi-cylinder 

n± = {(x,y) G Q : rkx > 0}. 

Take any A G C \ [0, oo) and consider the function 

(2.9) U = (L - A) _1 (F - AU) G H 2 (Q). 

The function W = U — Uq belongs to L 2 (Q) and satisfies the equation 
AW + AW = in each semi-cylinder Q±. Let W± stand for the restric- 
tion of W to Q,±. The functions W± can be expanded in the Fourier 
series 

(2.10) W±(x,y) =J2 w ne inv e~ lxlV ^ X , Re Vn 2 - A > 0. 

n 

Both series series converge in L 2 (Vt±) ) and 

/ \W±(x,y)\ 2 dxdy = V \ W p 2 

Jn ± ^2Re v / ^ r A 

Hence, W± G -^ 2 (^±) is equivalent to J2 n \w^\ 2 (n 2 + l) - ^ < oo. It 
follows that for each x G R the series in (2.10) converge in i/^^S 1 ) 
and moreover, W±(x, •) are continuous as functions of x with values in 
H-i^ 1 ). The same is true for the function U, and this explains the 
meaning of the condition (2.1): namely, 

(2.11) U(0+,y) = U(0-,y) as distributions in H'^(S r ). 

Denote by M(fl) the class of all functions U G L 2 (f2) which meet the 
following conditions. 

1. The distributions A(U \Q±) are functions in L 2 (Q±). 

2. The condition (2.11) is satisfied. 
For any A ^ [0, oo) we also set 

M A (Q) = [W G M(Q) : AW + AW = in Q ± }. 

The Fourier expansion of any function W G Ma(O) has the form 

(2.12) W(x, y) = J2 w n e my e^ xlV ^ X , 

n 
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that is, for the coefficients in (2.10) we have w+ = w n (= w n ) and thus 
the function W{x, •) is even in x. We also conclude that 

W G M A (Q) W(0, •) G if^S 1 ). 

Let us recall that in the terms of the Fourier coefficients w n the latter 
inclusion is equivalent to 

(2.13) ^ |w n | 2 (ra 2 + l) - * < oo. 

n 

Differentiation in (2.12) shows that for any W G Ma(^) the derivatives 
Wy(x, •), W' x (x, •) take values in the space H-2^ 1 ). The first of them, 
being an even function, is continuous in the topology of this space for 
all x G R 1 . The second one is continuous in the topology of /^^(S 1 ) 
for x > and for x < separately, and its jump across the circle 
{x = 0} is well defined as an element in /^"^(S 1 ). The decomposition 
U — Uq + W, where U is defined by (2.9), shows that the same is true 
for any U G In particular, this gives the precise meaning to 

both sides in (2.2) as distributions in //^(g 1 ). 

Substituting the Fourier expansion (2.12) and its differentiated forms 
into (2.2), we arrive at the system (2.4), (2.5), (2.6) which is equivalent 
to the initial problem. 

The following version of the Green formula is implied by the above 
argument. 

Lemma 2.1. For any U G M(Q) and V G H 2 (Q) (so that V(0,-) G 
if 2 (S 1 ) ) we have 



2. 




(AUV - UAV)dxdy 



= -[ (U x (0+,y)-U x (0-,y))V(0,y)dy, 

where the integrals on the left-hand side are understood in the sense of 
distributions on Q± and the integral on the right-hand side is understood 
in the sense of distributions on S 1 . 

Denote by 53 the differential operator appearing in the condition 
(2.2): 

(2.15) H>u — i(u' y cosy + (u cos y)' ). 

The operator 23 is symmetric as acting in the space L 2 (§ x ). The follow- 
ing useful equality, which is valid for U G M(f2) satisfying (2.2) with 
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an arbitrary a > and any V G H 2 (Q), is a direct consequence of 
Lemma 2.1: 



(2.16) 



(AUV - UAV)dxdy 



a / U(0,y)'BV(0,y)dy. 



s 1 



Indeed, substituting (2.2) into (2.14) and integrating by parts, we arrive 
at (2.16). 

3. The problem of self-adjointness 

In order to study self-adjoint realizations of L a for a > 0, we first 
of all introduce two sets, D Q and D*, on which the operator is well 
defined. It is convenient to do this in the terms of the expansion (2.3). 

Definition 3.1. An element U ~ {u n } G M(f2) lies in D a if and only 
if u n \R± G H 2 (R±) for all n e Z, the conditions (2.5) and (2.6) are 
satisfied, and 

/ \—u'n + n 2 u n \ 2 dx<oo. 

An element U E T) a belongs to D* , if the number of non-zero terms 
{u n } in the expansion of U is finite. 

We denote 

L a = -Af2) a , L* = -A CD*. 
Lemma 3.2. The operator L* is symmetric and 

= (KT- 

The proof is standard and we skip it. 

Theorem 3.3. 1) For < a < 1 the operator L a is self-adjoint and, 
hence, L* is essentially self-adjoint. 

2) For a > 1 the operator h a is non-self-adjoint, and the deficiency 
indices of~L' a are (2,2). 

Proof. We have to check whether the equation 

(3.1) L a W = {V a yW = AW 

with A ^ A has non-zero solutions W G D Q . If W is such a solution, 
then W G Ma(J1) and, by (2.12), each component in the expansion 
(2.3) for W can be written as w n (x) = w n e~\ x ^ n2 ~ A . The coefficients 
w n should satisfy conditions (2.6) that turn into 

(3.2) (n + l/2)w n+1 - 2a~ 1 w n Vn 2 - A + (n - l/2)iu n _i = 0. 
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The analysis of the system (3.2) is similar to the reasoning in [9], 
section 4, and is based upon the classical Birkhoff - Adams theorem, 
see [3], Theorem 8.36. The formulation of its leading case, which we 
need for the study of the operator ~L a with a ^ 1, is also reproduced in 
[9]. This theorem deals with one-sided sequences (n G N rather than 
n G Z as in our case), and we have to analyze the behaviour of w n for 
n — > +oo and for n — > — oo separately. 

For n — ■> +oo we find from the theorem that for a ^ 1 the equation 

(3.2) has two linearly independent solutions {w^(+)} such that 

(3.3) «;±(+) = (Aj) n n- 1 / 2 (l + 0(n- 1 )), A± = aT 1 ± Va- 2 - 1. 

For ra — > — oo we find in the same way that the system has two linearly 
independent solutions {w^(— )} such that 

(3-4) ' 

w ±(_) = (A ± )"|n|- 1/2 (1 + Odn]- 1 )), \t = -a' 1 ± Va~ 2 - 1. 

If a < 1, we conclude from the above asymptotic formulas that both 
for n > and for n < only one of the basic solutions decays as 
\n\ — > oo. Hence, the space of sequences {w n } satisfying (2.13) (or, 
equivalently, such that W G Ma(H)) is no more than one-dimensional. 
Suppose that {w n } is such a sequence, and apply the following identity 
for solutions of recurrence equations of the type 

Qn+lCn+1 + PnCn + QnC n -l = 0, Tl G Z, 

with Q n real, 
(3.5) 

N 

^2 \C n \ 2 \m.P n = -Q N+1 lm(C N+1 C N ) - Q_ Ar Im(C_ iV _ 1 C_ Ar ). 

n=-N 

The proof is straightforward and we skip it; cf. (4.23) and (4.24) in [9]. 
Applying (3.5) to the equation (3.2), we obtain 

N 

2a~ 1 |u; n | 2 Im \J n 2 — A = (N + 1/2) Im(wjv+iw]v + w_at_iW_at). 

n=—N 

By (3.3), (3.4) the right-hand side vanishes as — > oo. Since for non- 
real A the sign of Im y/n 2 — A is negative if Im A > and positive if 
ImA < 0, we conclude that w n = for all n G Z. It follows that for 
en < 1 the operator L Q is self-adjoint. 

If a > 1, then |A±| = 1 and by (3.3), (3.4) any solution {w n } satisfies 
(2.13). This shows that for a > 1 the operator L a is non-self-adjoint 
and the deficiency indices of L* are (2, 2). 
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Now, let a = 1. Then the case (ci) of the Birkhoff - Adams theorem 
applies, and the equation (3.2) has two linearly independent solutions 
of the form 

~ n ±v/3 ^, n — > oo 

and similarly for n — > — oo. For any non-real A only one of such 
solutions may satisfy (2.13). Using again the identity (3.5), we con- 
clude that the equation (3.2) has no non-zero solutions satisfying (2.13). 
Hence, the operator Li is self-adjoint. □ 



4. Using quadratic forms. Spectrum for a < 1 

For small a the simplest way to study the spectrum of the operators 
L a is to use quadratic forms. Our argument here follows the same line 
as in [9]. However, again, as in section 2, we have to take into account 
that the sequence {u n } is two-sided. 

Integrating by parts in the expression for (L a U, U) over the semi- 
cylinders il±, we find for U G £>*: 

(L a U,U) = [ \VU\ 2 dxdy 
Jn 

- [ U' x (0-,y)U{^y)dy+ [ U' x (0+,y)U(^y)dy. 

Taking into account the condition (2.2), we obtain 

(L a U,U)- [ \VU\ 2 dxdy 
Jn 

= ia I (C/;(0, y) cos y + (U(0, y) cos y)' ) U(0, y)dy 

= ia f {U' y (0,y)U{^-U(0,y)U^)) C osydy 
Js 1 

= -2af lm(U^O,y)U(0,y))cosydy. 

In the representation (2.3) this turns into 

(4.1) l a [U] := (L a U, U) = 1 [C/] - ab[U] 
where 

(4.2) l [U] = I (Kl 2 + n 2 \u n \ 2 )dx, 

(4.3) b[U] = J](2n-l)Re(« n (0K_i(0)). 
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Completing the set Do in the metric lo[U] + \\U\W 2 ^, we obtain a set 
which we denote by D. On d the quadratic form 1 is well defined and 
closed, and the associated self-adjoint operator in L 2 (Q) is L . Along 
with t), we need its subspace of co-dimension one, 

t)' ={U ~{t ln }6 0:n„(0) = 0}. 

Lemma 4.1. For any U G t)' the following inequality is satisfied: 

(4-4) \b[U]\<l [U]-\\u \\l 2{R) . 

Proof. Denote by c) + (by 0~) the subspace in D, formed by the elements 
U ~ {u n } whose all components with n < (with n > 0) are zeroes. 
For U e d ± we have b[C/] = ±b ± [f/] where 

(4.5) b^U] = J2( 2n - !) Re («±n(0)«±(„-i)(0)) • 

n>l 

The estimates for b + [U] and for b _ [?7] are identical and we carry them 
out for the 'plus' sign. We derive from (4.5) that 

\b+[U]\ < 1/2)(M0)| 2 + |^i(0)| 2 ) < ^2nK(0)| 2 . 

n>\ n>l 

Now to the n-th term in the last sum we apply the elementary inequality 

27l/(0)| 2 < / (|f| 2 + 7 2 m 2 )^, VfeH\R), 7 >0, 
Jr 

with 7 = n. We obtain 

\h + [U]\<J2 [ (\u' n \ 2 + n 2 \u n \ 2 )dx. 

n>l Jm 

Together with the similar inequality for b _ [C/], this yields (4.4). □ 

It is not difficult to show that the factor 1 in front of lo[U] on the 
right-hand side of (4.4) cannot be improved. 

With Lemma 4.1 at our disposal, it is easy to characterize the spec- 
tral properties of the operator L a for a < 1. 

Theorem 4.2. Let < a < 1. Then 

1) cr ess (L a ) = tr(L ) = [0, oo). 

2) The negative spectrum of~L a consists of exactly one non- degenerate 
eigenvalue. 

If a <\, then also 

3) <T a . c .(L a ) = <7 a . c .(L ) = [0,oo), m a . c .(L a ) = m a . c .(L ) 

(of. (2.8);. 
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The proof of the statements 1) and 3) basically repeats the argument 
in [9], section 9, and we skip it. To justify the statement 2), we first of 
all note that by Lemma 4.1, for a < 1 the quadratic form \ a , restricted 
to the domain D', is positive definite and closed. Since dimd/D' = 1, 
the quadratic form \ a , considered on the whole of D, is bounded from 
below and also closed. The corresponding self-adjoint operator is L a . 
For a — 1, the quadratic form l a is only closable on D, and the operator 
Li corresponds to the closure of li. 

This reasoning shows that for < a < 1 the number of negative 
eigenvalues of L a is no more than one. In order to show that it is 
exactly one, it is enough to find an element U E D which is such that 
la[U] < 0. To this end, we take U ~ {u n } with only two non-zero 
components uo,ui, then the desired inequality is 

Ool 2 + K| 2 + \ux\ 2 )dx < aRe(«i(0)ito(0)). 

It is satisfied, for instance, if we take U\{x) = e - ' 21 ' and u (x) = 
e~ x l 2 e~ e \ x \ with e sufficiently small. 

Remark 4.3. For a > 1 the quadratic form a Q is unbounded from below. 
We have to show that for any a > 1 and any M > there exists an 
element U ED, such that 

(4.6) aa^+MH^n^O. 

Choose a number N E N and take U ~ {u n }, where un(x) = e ~\ x W N ' 2 + M 

un-i(x) = e~\ x ^( N ~ lS>2+M and all the other components u n in the ex- 
pansion (2.3) are zeroes. Then 

/ (K| 2 + (N 2 + M)\u N \ 2 ) dx = 2VN 2 + M, 

[ (k'iv-il 2 + ((N - l) 2 + M)|m^| 2 ) dx = 2a/ (N — l) 2 + M, 

and 

a a [U] + M\\U\\ 2 L , {n) = 2(VN 2 + M+ y/(N - l) 2 + M) - a(2N - 1). 

It is clear, that for any a > 1 the last expression is negative, provided 
that N is taken large enough, and we are done. 

5. The case a > 1. Singular solutions 

In order to reach a better understanding of self-adjoint realizations 
of the operator L a for a > 1, we describe here the behaviour of the 
singular solutions found in section 3. 
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For a > 1 the asymptotic expressions for w„(±) as in (3.3) and (3.4) 
can be re-written in a simplified form. Indeed, set 

y(a) = arccos a -1 , 

then 

A+ = -Al = e iy(a \ Al = -A; = e' iy{a \ 

Therefore, 

u,±(+) = e ±iny ^n-^(l + 0(n' r )), n -> +oo; 



(5.1) 



«,±(-) = (-l) n e Tiny{a) \n\-^(l + 0(\n\- 1 )), v - - x . 



By (2.3), (2.12), and (5.1), each L 2 -solution of the equation (3.1) can 
be represented as 

W(x, y) = W(x, y; +) + + W(x, y; -), 

where K is a constant, W(x, y\ +) is a certain linear combination of 
the functions 



W ± (x,y;+) =^2w±(+)n-h iny e-W V ^ X 



n>0 



(5 2) 

= e ™(w±w(°)) n -§ e ^*l^/^(l + 0{n~ 1 )) , 

n>0 

and W(x, y; — ) is a linear combination of the functions 

W ± {x,y ] -)=Y,<< 
(5.3) " <0 



= ^ e m (f±fH— )| n |-5 e -l-lv / ^A(i + o(|n|- 1 )). 

n<0 

Note that y/n 2 — A = |n| + 0(n _1 ), and hence 

e -Nv^A = e -NH (1 + I^IOdn]- 1 )). 

Denote by V ± (x, y; ±) the functions obtained by replacing the factors 
e -\x\Vn 2 ~A e -|x||n| j n gg^h term of the sums in (5.2) and (5.3) and 
dropping the terms 0(|n| _1 ). The error is a bounded function rapidly 
decaying as \x\ — > oo. We have 

V ± (x,y;+) = J2n~h' n(lxhiiy±yia))) . 

n>0 

The behaviour of such sums as \x\ — i(y ± y{a)) — ► is well known. 
Say, it can be easily derived from the equations (13.11) in Chapter II 
of the book [10]. Denote 

z± + = \ x \ - i(y ± y(a)), z±~ = \x\ - i(y ± y(a) - tt), 
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then 

V ± (x, y;+) = C (z±+p + 0(1), z ± + - 0, 

with an appropriate choice of the branch of the square root, and some 
constant C. In the same way, 

V ± (x, y;-) = C (z ± -y* + 0(1), z±~ - 0. 

The reasoning above gives the following description of singular solu- 
tions W(x, y). These solutions depend also on the choice of A, but the 
leading terms of their singularities do not. For this reason we do not 
reflect dependence on A in our notations. 

Proposition 5.1. The singular solutions W(x,y) of the equation (3.1) 
have singularities at the points (0,yj) G Q, were yj, j = 1,2,3,4, are 
the points ±y(a) and ±y(a) + 7r(mod 2n). The singularity at each 
point is of the form 

W(x,y)~C j (\x\-iy j )--*+0(l). 

In order to explain the role of these four singular points, let us check 
the Shapiro - Lopatinsky criterion for the ellipticity of the boundary- 
value problem —AU = F under the conditions (2.1) and (2.2). In our 
case this criterion determines the point y G S 1 as regular if and only if 
the problem 

-<f>'\t) + 0(0 = 0, t^O, 0'(O+) -<j>'(0-) =±2acos#(0) 

has only trivial bounded continuous solutions on the line t G (— oo, oo). 
This requirement is violated exactly at the points y = yj, j = 1, 2, 3, 4, 
where a|cosy| = 1, the solution being <p(t) = e~l*L On the other 
hand, for a G [0, 1) the Shapiro - Lopatinsky condition is satisfied at 
all transition points. Therefore every weak solution of the equation 
— AU = AU satisfying (2.2) belongs to H 2 in both half-cylinders fl±, 
so it is non-singular, which explains the self-adjointness. 

6. The case a > 1. Spectral properties 

For a > 1, the main technical difficulty stems from the fact that 
Definition 3.1 does not describe the class T) a in the terms of standard 
function spaces on Q. For this reason, our argument here is rather 
lengthy. 

Let us fix some self-adjoint extension L a of the operator L*. The 
spectral properties discussed in this section do not depend on the choice 
of the extension. 
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We start by establishing a formula for the difference of resolvents of 
the operators L Q and L . The method for finding this kind of expres- 
sions is widely used and was proposed by Birman in [2]. Let first A be 
a non-real number. It belongs to the resolvent sets of both operators 
L Q and L , and we denote by R a , Ro the corresponding resolvents. 

Take some F,G G L 2 (il), and consider the sesqui-linear form 

(6.1) r[F, G] = ((R a - Ro)F, G) = (R a F, G) - (F, R* G). 
Denote 

R a F = U, RqG = V, 

then U G D a and V G H 2 (Q). Thus the quadratic form (6.1) can be 
re-written as 

(U,(L -A)V)-((t a -A)U,V)= ( [ +[ )(AUV-UAV)dxdy. 

\Jn + Jn-/ 

Applying (2.16), we arrive at 

r[F,G]=a [ U(0,y)'BV(0,y)dy, 

where B is the operator (2.15). Hence, the latter equality gives the 
representation of the operator R a — R as 

(6.2) R a - R = 2aS*T, T = TR a , S = SrR^, 

where T stands for the operator of restriction of functions on f2 to the 
circle x = 0. The operator T is bounded from L 2 (Q) to H~2 (S 1 ), and S 
is bounded from L 2 (Vt) to H 2 (S 1 ), so that S* is bounded from (S 1 ) 
toL 2 (Q). 

Our next step is to derive a pseudo-differential equation for the dis- 
tribution w = TW, where 

(6.3) W = U - V x := R Q F - R F, F G L 2 {VL). 

Evidently, W G Ma(Q) and thus w G H'^CB 1 ). Below we denote by A 
the operator —^2 in L 2 (S 1 ), extended to distributions on S 1 . It follows 
from the representation (2.12) that 

W^0+,y) - W' x $-,y) = -2j2^nV^Ae m y = -2(A - A)*w(y). 

n 

Now, taking into account the transmission conditions for U and for V±, 
we find that 



W^0+,y) - W x (0-,y) - aVWfry) = a'BV^y), 
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or 



(6.4) (2(A - A) 5 + qS ) w = -«SrR F G #5 (S 1 ) . 

The operator R a — R is, of course, bounded. We are going to show 
that, actually, it is compact. The proof is based upon the fact that 
the operator T in (6.2) acts from L 2 (Q) not only into if - 2(g) 1 ) but 
into a smaller space, iJ^S 1 ), for any e > 0. To show this, we need 
an a priori estimate for the equation (6.4). This equation is elliptic 
for a < 1, but for a > 1, which is the case we are dealing with, it is 
degenerate, so some more effort is needed. 

Lemma 6.1. For any e > there exist constants C,C such that for 
any w G H~2 (S 1 ) 

(6.5) \\w\\ H -^ } < C\\2(A - A)*w + aSHI^ipx) + C'\\w\\ H _ k§iy 
provided that the first term on the right-hand side of (6.5) is finite. 
Proof. Denote by P± the Riesz projections, 

p +f = ^ etky y ky > p -f= 77-1 E(/> eiky y ky - 

k>0 k<0 

Here the sums are understood in the sense of distributions; in particu- 
lar, if / G iP(S) 1 ), s G R, both series converge in iJ^S 1 ). 

The operators P± differ by smoothing operators from pseudodiffer- 
ential operators on the circle with symbols 

p+(y,v) = \] i 1 ! 7?> n' p-(y,v) = i-p+(y, , n), 

10 if i] < 0; 

see the discussion in [1] about the Fourier series representation of pseu- 
dodifferential operators on the circle. 

For w G if^S 1 ) we denote by w± the distributions w± = P±w. The 
operator (A — A) 2 is, up to a smoothing term, the pseudodifferential 
operator with symbol (rj 2 — A)2 = \r]\ + 0(\ri\~ l ). As it follows from 
the composition formulas for pseudodifferential operators in dimension 
one, the operators in (6.5) commute or almost commute with P±: 

(A - A)5 P ± = P ± (A - A)3, 'BP ± =P ± , B + K, 

with K being a smoothing operator. Thus, up to en error being an 
operator of order —1, the operator (A — A) 2 acts on the components 
w± as the differentiation, with proper coefficients: 

\\(A - A)^w ± T i<\\ H i m < C||^ ± ||^_i (sl) . 
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Therefore, (6.5) will follow as soon as we prove that 

(6.6) \\w ± U-^) < C\\ ±w> ± + ^^ ± ||^ (gl) + C'\\w ± \\ H _ hsl) . 

The estimate (6.6), even with — e replaced by 3/2 on the left-hand 
side, would follow automatically, if the operators ±2id y + aH> were 
elliptic for both signs ±. This is the case for \a\ < 1. However for 
| a | > 1 these operators have points of degeneracy of ellipticity, i.e. 
the points where the principal symbols (±1 + acosy)i] vanish. Note 
that these are exactly the points where the singularities of the singular 
solutions are located, see section 5. For such degenerate operators 
considering the principal symbol is not sufficient for getting a priori 
estimates, so the influence of lower order terms in £ must be taken 
into account. 

We concentrate on the case of the 'minus' sign in (6.6). Let us denote 
h{y) = a cosy and set 

, a , 1 

u = —w'_ + —H>w~ = (h(y) — l)w'_ + -h'(y)w-. 

We also set g = (h(y) — with a properly chosen branch of the 

square root. Note that g' = {h(y) — l)~5-u. Our next task is to derive 
an estimate of g in the terms of u, assuming that u G H^^Ei 1 ). 

The latter assumption on u implies that the function (h(y) — l)~^u 
belongs to the space H~ S (E> 1 ) for an arbitrarily small 5 > 0, say 5 < 1/2. 
To justify the above statement, we must show that 
(6.7) 



(%) - 1) 5 m(|/)C(|/)^ 



<^MI„V) IICII ^ S1 )' v Ce^(s x ). 



But this follows from the Holder inequality, since \h-l\-2 E //(S 1 ) for 
any r < 2, and by the embedding theorem u G /^(S 1 ) for any q < oo 

and C G L^S 1 ). 

It follows from (6.7) that 

\W\\h-H^) = ll(%) - l)"*u||tf-«(si) < C W u W H i (s iy 

Therefore, the function g = (h(y) — l)^W- lies in i7 1_<5 (§ 1 ) and satisfies 
the estimate 

ll5 r l| J F/ 1 -< 5 (s 1 ) < C W u W H h(si) + C "llfi'll^- JV (s 1 )' 
with iV being arbitrarily large. 
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By the definition of g, we have u>_ = (h(y) — 1) _5 <?. An estimate, 
similar to (6.7) (even a simpler one, since g G L 00 ^ 1 )), shows that u>_ 
belongs to H~ e ('S 1 ), with the required estimate. □ 

The estimate, just proved, enables us to establish the compactness 
of the difference of resolvents R a — R and of several related operators 
and prove spectral estimates. 

Proposition 6.2. The operator R a — R is compact, moreover for its 
singular numbers s n (R a — R ) the estimate 

(6.8) s n (R a - R ) = 0(n-^) 

holds for any e > 0. Further on, 
(6.9) 

s n ((R a - R )Ro) = 0(n-t +e ), s„(R (R a - R )) = 0{n~^). 
Proof. It follows from the factorization (6.2) that 

s n (R a -R ) <Cs n (T), 

where we have to consider the operator T as acting from L 2 (Vi) to 
H-h{^ 1 ). For F G L 2 {Vt) we define the function W as in (6.3) and take 
w = W{0,-), then 

TF = TR a F = w + TK F. 

The operator rR acts from L 2 (f2) to ifi(S 1 ), and the distribution w 
satisfies the equation (6.4), whose right-hand side belongs to if ^(S 1 ). 
Lemma 6.1 applies and gives w G if _e (§ 1 ). It follows that the operator 
T is bounded as acting from L 2 (Q) to //^(S 1 ), and therefore, the 
singular numbers of the operator T : L 2 (Q) — > H^2(E> 1 ) are controlled 
by those of the embedding if^S 1 ) -> iJ^^S 1 ). The latter are of the 
order 0(n~z +t ), whence the required estimate (6.8). 

Further on, we factorize the operator R (R a — R ) as 

Ro(Ra — Ro) = 2aRo5'*T '. 

Since we already know the singular numbers estimate for the operator 
T : L 2 {Vt) — > H^^i^ 1 ), it is sufficient for us to consider the operator 
RoS 1 * as acting between the spaces //^(S 1 ) and L 2 (Vl). It is more 
convenient to deal with the adjoint operator 

SR* = sr(R*) 2 : L 2 (n) -> ^(S 1 ). 

This operator is bounded as acting from L 2 {Vt) to if^S 1 ). Hence, the 
singular numbers of the same operator but considered as acting between 
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the spaces L 2 {Vt) and if 2 (S 1 ) are controlled by those of the embedding 
operator iff (S 1 ) -> ii^S 1 ). The latter are of the order 0(n~ 2 ). This, 
together with the estimate for T, proves the second estimate in (6.9). 
The first estimate in (6.9) follows from the second one by passing to 
adjoint operators. □ 

Now we arrive at our main result on the spectrum of the operator 
L Q , a > 1. 

Theorem 6.3. For a > 1 the spectrum of the operator L Q consists of 
the essential spectrum filling the semi-axis A > and the eigenvalues 
below the point 0. The set of eigenvalues below the essential spectrum 
is unbounded from below, may have onlyO and —00 as limit points, and 
for the counting function n(t) = #{A G o- disc (L a ), A G (—t, —to)}, with 
any fixed t > 0, the estimate holds 

(6.10) n(t) = 0(t 2+tl ), for any e l > 0. 

The absolutely continuous spectrum of L a fills the half-line A > and 
its multiplicity function coincides with that of L . 

Remark 6.4. The estimate (6.10) is rather rough. The authors believe 
that a more detailed analysis, based upon a further study of the de- 
generate equation (6.4), would show that the counting function has 
the asymptotics n(t) ~ Ctz as t — > 00. Moreover, we think that the 
negative eigenvalues do not have as their limit point. 

Proof. First, we note that due to Weyl theorem, the essential spec- 
trum of the operators R Q and Ro is the same, therefore the essential 
spectrum of L a coincides with that of L , so it is the half-line [0, 00). 
Thus, the spectrum of L a below may only consist of eigenvalues with 
possible accumulation points only at and —00. The latter point must 
be an accumulation point for eigenvalues since the operator L a is not 
semi-bounded from below, see Remark 4.3. The discreteness of the 
negative spectrum implies that there are real regular points of the op- 
erator L Q , these are all points below 0, which are not eigenvalues. We 
fix such regular A < and consider the resolvents R a , Ro at this point. 
Then the above construction of the operator R a — R and the estimate 
for its eigenvalues can be repeated, this time for the chosen real A. The 
spectrum of Ro coincides with the interval [0, —A -1 ], and 

R« = Ro + (Rq — Ro)- 

The operator R is non-negative, therefore, for any /i < the number 
of eigenvalues of R Q (counting multiplicities) in (— 00,/i) is not greater 
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than the number of eigenvalues of R a — R in the same interval. The 
latter quantity is estimated by means of the eigenvalue bound (6.8), 
which under an appropriate choice of e = e(ei) leads to (6.10), with 
t = -A and t = -(u _1 + A). 

In order to justify the statement on the absolute continuous spec- 
trum, let us consider the difference R„ — Rfj. We have 

Rq — R-0 = (j^a — Ro) 3 + R«Ro(Ra _ Ro) 

+R a (R a — Ro)Ro + Ro(R« — Ro) 2 + Ro(R-« ~~ R-o)) 

and due to the estimates (6.8) and (6.9) each term is trace class. By 
Kato's theorem, the absolute continuous parts of operators L a and L 
are unitary equivalent. □ 



7. An alternative model 

Here we briefly describe an alternative model, where a slight change 
in the setting leads to some major changes in the spectral behaviour. 
The family M a of differential operators acts on the strip Q' — R x (0, 7r) 
and is generated by the Laplacian —AU = —U" 2 — U y2 , the Dirichlet 
condition U(x, 0) = U(x,ir) = 0, and two additional conditions at 
x = 0: 

U(0+,y) = U(0-,y) (= 17(0, y)), 
U' x (p+,y)-U' x (p-y) = -ia (U' y (0, y) siny + (U(0, y) smy)' y ) , 
cf. (2.1), (2.2). 

The Fourier expansion for this case has the form 

^ /2 
U = 2_^u n (x)<f n (y), ip n (y) = d-smny 

n=l 

(in short, U ~ {u n }). The equation and the boundary and transmission 
conditions reduce to an infinite system of ordinary differential operators 
on the real axis, coupled by the conditions at x = 0: 

-AU ~ {-< + n 2 w n }, nGN; 

each u n is continuous at x = 0; 

<(0+) - <(0-) = ia((n + 1/2K +1 (0) - (n - l/2)u n ^(0)), 

with uq taken to be identically zero. 

For a = the system decouples, and we get an analogue of (2.7), 
but this time with summation over nGN. From here we derive that 
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the spectrum er(Mo) is absolutely continuous, fills the half-line [1, oo), 
and its multiplicity function is given by 

(7.1) m a . c .(A;M ) = 2[A], VA > 1. 

It is these two differences with L , the sequence of u n being one-sided 
and the spectrum of the unperturbed problem starting at 1 rather than 
at 0, that lead to the changes in the spectral properties of the perturbed 
operator. 

The study of the self-adjointness of M Q for a > follows the same 
line as for the operators L a in section 3. It turns out that the operator 
M Q , considered on the natural domain (cf. Definition 3.1), is self- 
adjoint for a < 1. If a > 1, the operator has a one-parameter family 
M a of self-adjoint realizations. The singular solutions, which define 
these realizations by v. Neumann's scheme, have two singular points 
(0, y^), with singularities of the order C(\x\ +i{y — y ± )y^ . The points 
y ± are the solutions of the equation asiny = 1, these are exactly the 
points where the Shapiro - Lopatinsky condition is violated. 

Similarly to the cylinder case, the spectral analysis of the operator 
M a for < a < 1 is based upon considering the quadratic forms. The 
quadratic form for M a is 

m a [U] := (M a U, U) = m [U] - ah[U] 

where 

™o[U] = J2 / {K\ 2 + n 2 \u n \ 2 )dx, 

neN ^ R 

b[U] = J2( 2n ~ l)Im(M0K^(0)), 

n>2 

cf. (4.1), (4.2) and (4.3). The quadratic form m is positive definite 
and closed on its natural domain which we again denote by 0. The 
associated self-adjoint operator in L 2 (Vt') is M . The inequality 

(7.2) \h[U]\ <m [U], Uet, 

is checked in the same way as (4.4), and this time no second term as 
in (4.4) appears. The constant factor 1 in the estimate (7.2) is sharp. 
Hence, for a < 1 the quadratic form m Q is positive definite and closed 
on D. The corresponding self-adjoint operator in L 2 (Q') is M a . It is not 
difficult to show that for a > 1 the quadratic form m a is unbounded 
from below. 

We pass now to the description of the spectrum of M Q . It is here 
where the differences with L Q manifest themselves, cf. Theorem 4.2. 
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Theorem 7.1. Let < a < 1. T/ien 

j; <7 ess (M Q ) =<t(M ) = [l,oo). 

(J a . c .(M a ) = a a . c .(M ) = [1, oo), m a . c .(M Q ) = m a . c .(M ) 

(of- (7.i);. 

3) The spectrum o/M Q fre/ow; i/ie threshold Xq = 1 is finite. 

We skip the proof which basically repeats the argument in [9], section 
9. Note that one can also prove that for the pairs M Q , M and M , M a 
there exist complete isometric wave operators. 

The quadratic form mi \ d is non-negative and closable, it generates 
the operator Mi. It is possible to show that its essential spectrum is 
the half-line [0, oo). 

The analysis of the discrete spectrum of M a for a G (0, 1) is based 
upon a version of Birman-Schwinger principle found in [9]. Before giv- 
ing its formulation, let us recall the following well-known notations. 
Given a real number A and self-adjoint operator Q, whose spectrum on 
(— oo, A) is discrete, we write N_(X; Q) for the number of the eigenval- 
ues X n (Q) < A, counted according to their multiplicities. We also write 
N + (X;Q) = N.(-X;-Q). 

It turns out that within an error which is no greater than 1, the 
number JV_(l;M a ) coincides with jV + (a: -1 ; J), where J is a certain 
infinite Jacobi matrix: 

(7.3) < JV_(l;M a ) -N+ia- 1 ; J) < 1. 

The reasoning is the same as in [9], however the Jacobi matrix J turns 
out to be different: it is the zero-diagonal Jacobi matrix, with the 
non-diagonal entries given by 

n_1 / 2 

2jn,n-l - 2jn-l,n ~ (^2 _ 1)1/4(^2 _ 2n )l/4 ' 

Since j n ,n-i 1/2 as n — > oo, the matrix J has the absolutely con- 
tinuous spectrum filling the segment [—1, 1] and the spectrum outside 
this segment is discrete. Note that a G (0, 1) is equivalent to a^ 1 > 1, 
so that both terms in (7.3) are finite. 

In order to estimate N + (/i;J), \i = a -1 , we use the asymptotics of 

Jn,n— 1 • 

(7.4) j n ,n-i ~ ^ + ^n~ 2 + o(n' 2 ), n -> oo. 

Using the results of Geronimo [6], [7], combined with some standard 
variational tools, one can show that jV + (//;J) can be estimated from 
below and from above by | log(/i — 1)|, with different constants. Thus 
the number of eigenvalues of M a in (0, 1) grows logarithmically as 
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a / 1. We believe that actually a logarithmical asymptotics for the 
eigenvalues holds. 

When a becomes larger than 1, the phase transition occurs, similar 
to the cylinder case. Each self-adjoint realization M a of the operator 
M a is unbounded from below, with the spectrum below the point 1 
being discrete. The absolutely continuous spectrum is still the half- 
line [l,oo), with the same multiplicity function as for M . All these 
properties are proved using the methods exposed in section 6. Some 
additional technical complications are caused by the fact that now we 
should prove estimates of the type (6.5) for the operators on an in- 
terval (0,7r), rather than on the circle S 1 which is a manifold without 
boundary. But these complications can be overcome. 
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